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Abstract. A version of the saddle point method is developed, which allows one to describe 
• exactly the asymptotic behavior of distribution densities of Levy driven stochastic integrals with 

\ deterministic kernels. Exact asymptotic behavior is established for (a) the transition probability 

Ph ' density of a real-valued Levy process; (b) the transition probability density and the invariant 

\ distribution density of a Levy driven Ornstein-Uhlenbeck process; (c) the distribution density of 

the fractional Levy motion. 

•)—> . 

I. Introduction 

In this paper, we develop a version of the saddle point method, which allows one to describe 
p\j exactly the asymptotic behavior of distribution densities of Levy processes and, more generally, 
^ \ Levy driven stochastic integrals with deterministic kernels. We start the exposition with the 
qq \ outline of the principal idea of the approach. 

Let (Z t ) t >o be a real-valued Levy process with characteristic exponent ip; that is, 

^ (1.1) Ee izZt = e^ iz \ t>0. 

The function if> : R — > C (the characteristic exponent of the process Z) admits the Levy-Khinchin 
representation 

•rH . „ 

X (1.2) if;(z) =iaz-bz 2 + / (e iuz - 1 - izul { \ u{ < 1} ) ^(du), 

g ■ Jm 

where a G R, b > 0, and //(•) is a Levy measure, i.e. f R (l Au 2 )fi(du) < oo. Under some conditions 
(see Section 2 below), the function e 1 ^ is integrable, and hence the transition probability density 
Pt{x) of the process Z t has the integral representation as the inverse Fourier transform of the 
characteristic function (1.1): 

(1.3) Pt(x) = ±- [ e~ izx+t ^dz. 

Our intent is to investigate the oscillatory integral (1.3) using the saddle point method. According 
to this method (see [Co65]), one can, under the assumption that the characteristic exponent if) 
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admits an analytic extension to the complex plane, apply the Cauchy theorem in order to change 
the integration path in (1.3): 



(1.4) 



Pt 



-izx+tip(z) ^ 



Here C is certain properly chosen contour that allows one to apply the Laplace method ([Co65], 
[Ev62], [Fe77]) for estimating integral (1.4). A perfect choice of the contour C would be the proper 
branch of the curve {z : Im(— izx + tip(z)) = Im(— izoX + tip(z ))}, where Zo is a critical point of the 
function — izx + tip(z) (a saddle point). Under such a choice the integrand in (1.4) is real-valued; 
in this case the saddle point method coincides with the fastest descent method, see [Co65]. However 
the complicated "oscillatory" structure of the Levy-Khinchin representation of ip does not give an 
opportunity to solve the equation Im(— izx + tip(z)) = Im(— iz x + tip(z )) explicitly. Instead, we 
put in (1.4) C — K. + i^Q with i£ being a critical point of the function —izx + ttp(z). Under such 
a choice, we develop an appropriate version of the Laplace method and give exact asymptotics for 
the transition probability density pt{x). 

The saddle point method is a classic tool for estimating a distribution density in various versions 
of the local limit theorem with the normal domain of attraction (see [IL65], chapters 8, 10, and 
the references therein). In the Levy processes setting, the idea of applying the complex analysis 
technique was used, for instance, in [KS08] for getting upper estimates for (1.3) in the case when 
the characteristic exponent is real valued. 

Since we require the characteristic exponent ip to have an analytic extension to the complex 
plane, a standing assumption on the Levy measure within our approach is that it is exponentially 
integrable; that is, 



Equivalently, (1.5) means that the variable Z 1 has exponential moments, i.e. Ee cZl < oo for all 
cGR, see [Sa99], §25 — 26. Assumption (1.5) is non- restrictive, and is satisfied, for instance, for 
a generalized tempered Levy measure of the form fi(du) = ip(u)jl(du), where fx is another Levy 
measure, and ip has a super-exponential decay, i.e., e Cu ip(u) — > 0, u — > oo, for all C G K. For 
various results on generalized tempered Levy processes and models that lead to processes of such 
a type, we refer the reader to Rosinski and Singlair [RS10], Sztonyk [StlOa], [StlOb], Bianchi et. 
al. [BRKF08], [KRBF09]. The notion of a generalized tempered Levy measure is closely related 
to the notions of a tempered and a layered Levy measure, with the function ip{u) in the above 
definition respectively being completely monotonous or having a polynomial decay rate (for both 
these classes (1.5) fails). For the results on tempered and layered Levy processes and related 
models, see Rosinski [Ro07], Cohen and Rosinski [CR07], Cont and Tankov [CT04], Carr et. al. 
[CGMY02], [CGMY03], Baeumer and Meerschaert [BM09], Kim et. al. [KRCB09], Houdre and 
Kawai [HK06]. Of course, this list of references is far from complete. 

The method described above can be extended naturally for Levy driven stochastic integrals 
with deterministic kernels. Let 



(1.5) 




e Cy n{dy) < oo for all Cel. 



(1.6) 
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where I C R is an interval, / is a deterministic function, and Z t is a Levy process (in some 
particularly important cases, one should take / = R, and then Z should be assumed to be two- 
sided; see details in Section 2 below). The characteristic exponent of Y t can be written explicitly 
(see (2.4) below), which makes it possible to apply the method described above to study the 
asymptotic behaviour of the distribution density of Y t . 

We mention two particular classes of processes, frequently used in applications, and having 
representation (1.6). The Levy driven Ornstein-Uhlenbeck process is defined as the solution to the 
linear SDE 

(1.7) dX t = 7 X 4 dt + dZ t , t > 0, 

and has the integral representation 



1.8) X t = e 7 % + I e 7( *" s) dZ s , t > 0. 

Jo 



If the initial value X is non-random, the distributional properties of X t are determined by the 
second term in the right hand side of (1.8), which clearly has the form (1.6) with I = R + and 
f(t,s) = e 7 ^ _s ^l s < t . In what follows, we call such a process a non- stationary version of the 
Ornstein- Uhlenbeck process. 

The Ornstein-Uhlenbeck process is Markov one. It is ergodic (i.e. possesses unique invariant 
distribution), if and only if, 7 < and 

(1.9) / In \u\fi(du) < +00; 

J\u\>l 

see [SY84]. Clearly, our standing assumption (1.5) provides (1.9). Respective stationary version 
of the Ornstein- Uhlenbeck process can be represented as 



x t = f e^~ s Uz s , te 

J —00 



which is clearly of the form (1.6) with I = R, f(t,s) = e 7 ^~ s - ) l s < t . Conditions on the existence 
and smoothness of the distribution densities for Levy driven Ornstein-Uhlenbeck processes were 
studied in [Ma04], [BK08], [PZ09], [S10]. In some exceptional stationary cases, the density can be 
represented explicitly, see [BNS01]. However, as far as we know, any references concerning general 
estimates or a description of the asymptotic behaviour of such a density are not available. 

Another example of a process of the type (1.6) is the fractional Levy motion, defined, analogously 
to the fractional Brownian motion, by the stochastic Weyl integral 

(1.10) z H (t) = T{H l 1/2) I [(t - sf + - 1/2 - {-s) H + - l/2 \ dz s , t e R, 

where x+ = max(x,0), and H G (0, 1) is the Hurst index; see [ST94], [BCI04], [M06], [KM10] and 
references therein. In what follows, we will study the asymptotic behaviour of the distribution 
density of Z H {t) under the assumption that H > 1/2, which is the so called long memory case, 
see Definition 1.1 in [M06]. Note that in this case Z H is not a Markov process, in contrast to the 
Levy process Z, or the Levy driven Ornstein-Uhlenbeck process (1.8). 
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Heat kernel estimates for symmetric jump processes were studied systematically by Barlow, 
Bass, Chen and Kassman [BBCK09], Chen, Kim, Kumagai [CKK08], [CKK10], [C09], Barlow, 
Grigoryan, Kumagai [BGK09], Chen, Kumagai [CK03], [CK08], Chen, Kim, Kumagai [CKK09]; 
see also Bass and Levin [BL02] for the transition density estimates for a Markov chain on Z d . 
The approach used in the papers listed above relies on the paper by Carlen, Kusuoka and Stroock 
[CKS87]. For heat kernel estimates in domains we refer to the papers by Bogdan and Jakubowski 
[BJ07], Banuelos and Bogdan [BB04], Bogdan, Grzywny, and Ryznar [BGR10], Chen, Kim and 
Song [CKSa] - [CKSd]. Of course, this list of references is far from complete. 

In particular, heat kernel estimates for symmetric jump processes on M. d with jump kernel 
J(x,y), either bounded both from above and below by ^ x _^ d+a l|g-y|<i; < a < 2, d > 1, 

or decaying as e~^ x ~ y ^ , (3 e [0, oo), as \x — y\ — > oo, are studied in [CKK08] and [CKK10], 
respectively. Under the particular choice of the jump kernel J(x,y) = J(x — y), the processes 
studied in [CKK08] and [CKK10] become symmetric Levy processes. We postpone to Section 3 
(Example 3.4) the detailed comparison of the asymptotic results for the distribution densities of 
such processes obtained in [CKK08] and [CKK10], with the results obtained by our approach. 
Here we just mention that our approach, based on the the complex analysis technique, can be 
applied both for non-symmetric Markov jump processes, like the Levy driven Ornstein-Uhlenbeck 
process, and for non-Markov processes such as the fractional Levy motion. 

Let us outline the rest of the paper. Our main result on asymptotic behavior of the distribution 
densities of Levy driven stochastic integrals Y, including the Levy process Z itself, is formulated 
and proved in Section 2. To simplify the exposition, we give one-sided asymptotics; that is, we 
formulate the main result for the distribution density pt(x) only for x > 0. Clearly, one can easily 
deduce from this result the two-sided asymptotics, assuming additionally that the Levy measure 
of the process (-Z) satisfies conditions of Theorem 2.1. 

Conditions of the main result, Theorem 2.1, are quite abstract, and require an additional analysis 
in order to provide a verifiable criteria. For the reader's convenience and to clarify the exposition, 
we separate such an analysis in two parts. In Section 3 we consider an "individual" asymptotic 
behavior of the distribution density of Y t with fixed t. We formulate an "individual" version of 
Theorem 2.1 with verifiable conditions on the Levy measure \x and the kernel /. These conditions, 
in particular, reveal the "smoothifying" effect provided by the kernel /: typically, both to provide 
existence of the distribution density of the Levy driven stochastic integral Y t and to describe its 
asymptotic behavior, fewer restrictions on the Levy measure are required than in the case of the 
Levy process Z t itself. An illustrative example of such an effect is provided by the fractional Levy 
motion, where the assumptions on the Levy measure are finally reduced to 



In Section 4 we establish the asymptotic behavior of the distribution density of Y t , involving 
both state space variable x and time variable t. To shorten the exposition, we restrict ourselves 
to the case of a self-similar kernel /. The class of the Levy driven stochastic integrals with 
self-similar kernels, although not being the most general possible, is wide enough to cover the 
important particular cases of the Levy process Z itself and the fractional Levy motion Z H . As a 



(1.11) 
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corollary of the main result of Section 4 (Theorem 4.1), we obtain asymptotic relation 



'1.12) p t (x) = e Wz \V, t + x^oo, (t,x) G [* ,+oo) x 

hntJCz (f ) 



for the distribution density of the Levy process Z, and 

(1.13) p t (x) 1 e w zH(jHfrj2) ? t + x ^ 00 ^ (t, x) G [to, +oo) x 

for the distribution density of the fractional Levy motion Z H . Here t > is arbitrary, V Y , JCy 
with Y — Z, Z H are some functions, defined in terms of the Levy measure /i and the kernel /; 
see Section 4 below. Observe that the asymptotic formulae for distribution densities of Z and 
Z H possess the self-similarity property in spite of the fact that, in general, the families of these 
densities are not self-similar. 

Formally, Z H includes Z as a partial case with H = 1/2, and (1.13) with H = 1/2 transforms 
to (1.12). However, there is a substantial difference between the conditions under which these 
asymptotic results are available (see Corollary 4.1). To get (1.12), one should impose some "regu- 
larity" conditions (iVi) and (C) together with some "tail" conditions (Ti) and (T2). To get (1.13) 
with H G (1/2,1), it is sufficient to claim only "tail" conditions and non-degeneracy condition 
(1.11): there is no need for additional "regularity" conditions. Such a difference is caused by the 
"smoothifying" effect provided by the kernel in the integral (1.10). 

Theorem 2.1 and Theorem 4.1 describe the asymptotic behaviour of the distribution density 
precisely, but in an implicit form. In Section 5 we use these theorems in order to deduce explicit, 
although less precise, asymptotic expressions. In the same section we give another application of 
Theorem 2.1, and study the asymptotic behavior (as x — > 00 for a fixed a) of the ratio 

1-14 r a {x) = Fy ' 

p(x) 

for the invariant distribution density p of the Ornstein-Uhlenbeck process. Such a study is of 
particular theoretical interest, since the ratio (1.14) appears in the formula for the generator of the 
dual (i.e., time-reversed) process corresponding to the solution to SDE (1.7). Therefore, knowledge 
of the asymptotic properties of (1.14) would be useful when one is interested in studying the 
stationary version of the solution, respective Dirichlet form etc. For instance, in the forthcoming 
paper [Ku09] the estimate given in Theorem 5.2 below is used substantially in the proof of the 
spectral gap property for the Levy driven Ornstein-Uhlenbeck process. 

Formula (1.12) and Theorem 5.1 provide a detailed description of the asymptotic behavior of the 
distribution densities of the Levy process and the fractional Levy motion. This behavior exhibits 
two different regimes. In the first regime, where the ratio | (resp., tH + 1/2 ) stays bounded, the 
principal behavior of pt(x) is determined by the values of the functions XV, Ky (with Y = Z or 
Z H ) on a bounded domain. For instance, for any x > 

(1.15) ft (fx) ~ -_I== e «>*M t ^ +OQ 
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(for the Levy process Z) and 

(for the fractional Levy motion Z H ). In the second regime, where the ratio | (resp., ^172 ) tends 
to +00, the principal behavior of pt(x) is determined by the asymptotics of XV, /Cy (with F = Z or 
on +00. Such asymptotics are described in Theorem 4.1 for two cases: for the Levy measure 
\x being either "truncated" (i.e. supported in a bounded set) or "exponentially damped" (i.e. its 
tail satisfies certain exponential estimate, see (3.23)). This description gives some constant c*, 
determined in terms of the Levy measure \x only (see (5.15) and (5.16)), such that the statements 
below hold true (see Corollary 5.1 and Corollary 5.2 below). 

I. Case of the Levy process Z . For any constants c\ > c* and c 2 < c* there exists y = y(ci, c 2 ) 
such that for x/t > y, either 

(1.17) exp ^— cirrln ^— j j < pt(x) < exp ^—c 2 a; ln 
(if fj, is truncated), or 

(1.18) exp ^-cixln^" (j-^ < Pt(x) < exp (-c 2 xlnV" (^j^j 

(if \i is is exponentially damped). 

II. Case of the fractional Levy motion Z H . For any constants c\ > c* and c 2 < c* there exists 
y = y(ci, c 2 ) such that for x/t H+1 / 2 > y, either 

(1.19) 

ex K"TWTW^ ln (^)) ^ ^ ^ exp (~ f(g + vV-V2 ln (f^)) 

(if /i is truncated), or 
(1.20) 

/ C\X P=± ( X \\ ( C2X £=± ( X \\ 

exp (,"r(^ + 1/2)^-1/2 ln ' J ^ ^ ^ exp (,"r(^ + 1/2)^-1/2 ln K^wi) ) 

(if /x is is exponentially damped). 

In this paper we restrict ourselves to the case of one-dimensional processes in order to make 
the exposition reasonably short, and to give the main results in their most transparent form. 
These results have straightforward generalizations to the multi-dimensional case; we postpone the 
discussion of these generalizations to a further publication. We also restrict our considerations of 
the Levy process Z and the fractional Levy motion Z H to the case where the time variable t is 
separated from 0. The small time estimates require additional analysis of the local behavior of 
the Levy measure of the noise; this analysis is performed in the separate article [KK11]. 

2. The main result 

2.1. Preliminaries. Everywhere below Z is a Levy process and ip is its characteristic exponent; 
that is, (1.1) and (1.2) hold. 
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To exclude from consideration the trivial cases, we assume that b = and /x(R) > 0; that is, 
Z does not contain a diffusion part, and contains a non-trivial jump part. Moreover, we assume 
that /i satisfies (1.11), which is motivated by our intent to analyze the distribution density on 
the positive half-line. Finally, we assume Z to be centered, which means that the characteristic 
exponent is of the form 

(2.1) ij;(z) = [ (e iuz - 1 - izu) n(du), zeR. 

Jr 

This assumption does not restrict the generality: under (1.5), the increments of Z have moments 
of all orders, therefore the difference between the processes with characteristic exponents (1.2) and 

(2.1) is given by the explicitly calculable constant, which clearly does not effect the distributional 
properties. 

We consider Levy driven stochastic integrals of the form 

(2.2) Y t = Jf(t,s)dZ s , teT, 

where Tclis some set, and I C R is an interval. We allow the case where the interval I belongs 
not only to the half-line, but to whole R. In this case, the process given by (2.2) is assumed to be 
well defined on the whole line R, and to have independent and stationary distributed increments, 
with the characteristic exponent of the increments still being of the form (2.1). A standard version 
of such a process is the so called two-sided Levy process 

Zl t > 
-Z 2 _ t _, t<0' 

where Z 1 and Z 2 are two independent copies of a Levy process, defined on R + . 

We interpret (2.2) as an integral with respect to an infinitely divisible random measure; for the 
general theory of such integrals we refer to [RR89]. The integral (2.2) is well defined if, and only 
if, 

(2.3) Jf{t,s)ds < +oo, teT, 

and in that case its characteristic function admits the representation 



(2.4) Ee izYt = exp 



[ [ ( e l2/(M)u - 1 - izf(t, s)u) fi(du)ds 
J i Jr 



z e R, teT, 



see Theorem 2.7 from [RR89]. In what follows we assume that / satisfies (2.3), and /(t, •) is 
bounded for every teT. To exclude the trivial case Y t = a.s., we assume f I f 2 (t,s)ds > 0, 
teT. We also assume 

(2.5) J(f(t,s)V0) 2 ds>0, teT. 

This does not restrict generality since otherwise one can consider —Y t instead of Y t . 
For a Borel set icK, denote 

Q(t,z,A)= (1 - cos( f(t,s)zu))/i(du)ds, teR + ,^eR. 

J J{(s,u)£lxR: f(t,s)u<=A} 
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The functions G(-, -,A), with properly chosen sets A, will be used below as a tool for studying 
the properties of distribution densities of Levy driven stochastic integrals Y . One statement of 
such a type is formulated in the proposition below, which is in fact the classic Hartman-Wintner 
sufficient condition ([HW42]), reformulated in the context of Levy driven stochastic integrals. 

As usual, we denote by C^{R) the class of function, continuous and bounded together with their 
derivatives up to order k. 

Proposition 2.1. For given t G T, k G Z +; and \z\ large enough, let 

(2.6) B(t,z,R) > (k + l + 5)\n\z\ 
with some 5 > 0. 

Then Y t has a distribution density p t , which belongs to the class C*(R). 
In particular, if for a given t G T 

(2.7) &(t,z,R) ^>\n\z\ as \z\ ->■ oo, 
then Y t has a distribution density p t G C£°. 

Proof. By (2.4), condition (2.6) implies 

(2.8) \Ee izYt \ < \z\ ~ k ~ 1 ~ 5 for \z\ large enough. 

Hence the required statement follows by the inversion formula for the Fourier transform. □ 

As usual, we write /(£) ~ g(£), £ ->■ oo, or /(£) = o(#(£)), f ->• oo, if lim^^ ^ = 1 
or lim^oo 4|y = 0, respectively. We also use the notation /(£) <?(£)> C — ^ instead of 

/(0 — (9(0)i £ — ^ 00 ) when it is more convenient. The same conventions are used when 
functions / and g depend on t and/or on x. 

2.2. The main result: formulation and discussion. Since / is bounded, the exponential 
integrability assumption (1.5) implies that for t G T the function 

(t, z) = I I (e~ izf{t ' s)u - 1 + izf{t, s)u) fi(du)ds, t G T, z G C, 
J 1 Jr 

is well defined and analytic with respect to z. Denote 

H(t, x, z) = ixz + ty(t, z), 
and observe that, assuming (2.6), we have 

(2.9) p t (x) = — [ e H(t ^ z) dz, xeR, 

2tt J r 

which is just the inversion formula for for the characteristic function of Y t , combined with the 
change of variables z 1— > —z. 

Denote 

M k (t^) = —^(t,i£), k>l, ieR. 
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Clearly, 

M k (t,Z)= f [ u k f\t,s)e^ u v(du)ds = -^H(t,x,tO, k>2. 
J i Jm. ai i 

Since fi and f(t, •) are assumed to be non-degenerate, we have A^t^O > 0. Therefore there 
exists at most one solution £(t,x) to the equation 

(2.10) ^H(t,x,i£) = 0. 
Clearly, for any t G T we have £(£, 0) = 0. Note that 

■Mi(f , = £z*(t, x, i£)= f j uf(t, s) (e^ u - l) iM(du)ds = [ - l)ihj(dv), 

where /i t j denotes the image of the measure fi(du)ds under the mapping 

/xK9(s,m)4 f(t,s)u G K. 

Under the assumptions (1.11) and (2.5), which we assume to hold everywhere below, we have 
/j, t j(R + ) > 0. Therefore Mi(t,£) — > +oo as £ — > +oo, which means that is well defined 

and positive for x > 0, and 

(2.11) x) — >■ +00, as x — > +oo. 

Note that z = i^(t,x) is the unique critical point for H(t,x, •) on the line iR. 
We put 



^2 

D{t, x) = H(t, x, i£(t, x)), K(t, x) = M 2 (t, £(t, x)) = g^H(t, x, i£) 



In the sequel, we fix A C T x R + and denote 

T={t:3xe R + , (t, x) G A}, B = {(t, : 3(t, x) G .A, (t, £) = (i, a;))}. 

For instance, if .4. = T x M+ with some T C T, then T = T' and B = A. 

In the following theorem, which represents the main result of the paper, the function 6 : T —f 
(0, +oo) is assumed to be bounded away from zero on T, and the function x '■ T — > (0, +oo) is 
assumed to be bounded away from zero on every set {t : 6{t) < c}, c > 0. For a particular process 
Y, the choice of the "scaling" functions 9 and x is determined by the structure of the kernel /, 
see Section 4 below. 

Theorem 2.1. Assume that the following conditions hold true: 

(H x ) M 4 (t,C)<:M 2 2 (t,0, 0{t) +^oo, (t,0eB. 

K(^^tl) vi ) +in (( in ^ o ) vi ) 

«ln0(t) + x(t)£, ^(t) + £^oo, (U)eB. 

(H 3y ) T/iere exist R > and 5 > snc/i £/ia£ 
(2.12) 6(M,R + ) > (1 + 6)ki\z\, teT, \z\>R. 
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2 A1). 



(H4) There exists r > such that for every e > 0, 

inf Q(t,z,[r X (t), +00)) >6(t)((e X (t)) 

teT,\z\>e \ 

Then for every t e T £/ie /aw o/F t /ias a continuous bounded distribution density pt(x), and 
(2.13) p t (x) ~ ^_L= e D ^>, 9(t)+x^ 00, (t,x)eA. 

Remark 2.1. (On conditions). The conditions of Theorem 2.1 are rather technical and abstract. 
In Sections 3 and 4 below we give their more explicit versions, formulated in terms of the Levy 
measure fi and the kernel /. Note that (Hi) and (H 2 ) are, in fact, the assumptions on the growth 
of the tails of the Levy measure [i. In addition, (H2) is balanced with (H4), which in turn is 
closely related to the so called Cramer condition (see, for example, [Lu79], [IL65], Chapter 3 §3, 
and the discussion prior to Lemma 3.5 below). Finally, (H 3 ) is a proper uniform version of the 
Hartman-Wintner condition, see Proposition 2.1. Clearly, one can consider the stronger version 
of condition (2.12) with k + 1 + 8 instead of 1 + 5 (cf. (2.6) and (2.7)), and provide the asymptotic 
relations similar to (2.13) for the derivatives of the distribution density pt(x) up to order k. 

Remark 2.2. (On relation (2.13)). 1. Note that the asymptotic relation (2.13) corresponds 
completely to the standard form of an asymptotic relation obtained by the Laplace method. 
Typically, within this method one can prove that the integral 



/ 

J (a, 



e- F{x ^ dx 



is asymptotically equivalent to 



/ 2-7T 

(2.14) -e~ F ^\ x A :=argminF(A,x). 

V ^xx\^ X X) 

Clearly, (2.13) is exactly of the form (2.14) with appropriate F and additional normalizer l/(27r), 
which comes from the inverse Fourier transform formula. 

2. Our approach is in some sense related to the Large Deviations Principle (LDP). Namely, if 
Pl(dx) is the probability measure associated with Y t l := y Yl\=\ %h where {Zl} l i=1 are independent 
copies of (Z t ) t >o, then P}(dx) satisfies the LDP with a good rate function A t (x) := —D(t,x), in 
the sense that for all measurable subsets A C R 

1 1 
inf A t (x) < lim inf - In P\ (A) < lim sup - In P\ ( A) < - inf A t (x) ; 

xGinterior(A) I— >oo / i^oo ' a;Gclosure(A) 

see [DS89]. Moreover, assuming the exponential integrability condition (1.5) and existence of the 
transition probability density pt(x) for t > to, it is shown in [KS08] that 

(2.15) \im lnp f x) =D(t,x), 
cf. (1.15) and (1.16). 



10 



2.3. Proof. Note that Q(t,z,A) depends on the set A monotonously. Hence (#3) yields (2.6), 
and therefore (2.9) holds. In what follows, we analyze the right hand side of (2.9). We divide this 
analysis into several steps. 

Step 1: changing the integration contour. We prove that 

(2.16) p t {x) = — [ e H ^dz = — [ e H(t*,v+it(trt) dr] 

Recall that we assumed x > 0, which in turn implies £(t,x) > 0. Consider the domain 

(2.17) G M :={zGC:ImzG[0,((t,4 Rez 6 [-M, M], M > oj. 
The function H(t,x,z) is analytic in Gm, hence by the Cauchy theorem 

(2.18) / e H(t ' x ' z) dz = 0. 

JdG M 

Consider the integrals 

(2.19) f e H(t,x,±M+iv£(t, X )) dv _ 

Jo 

We have 

ReH(t,x,r] + iZ) = -x£- [ [ (l - e /(M)?u cos(/(t, s)rju) + f(t, s)£u)fi(du)ds 

J 1 Jr 

(2.20) =H(t,x,i£)- f [ e f(t ' s ^ u (l -cos(f(t,s)rju))fi(du)ds, ^el. 

J 1 Jr 

The function £ h-> H(t,x,i£) is real- valued, convex, and attains its minimal value at the point 
£(t,x). Then H(t,x,i£) < H(t,x,0) = for £ e [0,£(t,x)}. On the other hand, for every £ > 

I [ e nt ^ u (l - cos(/(f , s)r)u))v(du)ds > ff (1 - cos(/(t, s)r)u))n{du)ds 

J 1 Jr J J {(s,ti)e/xi: /(t,s)u>o} 

= e(t,?7,R + ). 

Therefore 

Re#(t,a;,±M + iuf(*,a;)) < -Q{t,±M,R + ), v e [0,1]. 

Thus, condition (H 3 ) implies that the integrals in (2.19) tend to as M — > +00, which together 
with (2.17) gives (2.16). 

In what follows we denote 

R(t, x, rf) = Re H(t, x, 77 + i£(t, x)), 

I(t, x, 77) = lmH(t, x, 77 + i£(t, x)) — xr] — t / {e^ l ' x ^ u sin(rju) — r/ujn(du). 

Jr 

Since a distribution density is real valued, we derive from (2.16) 

(2.21) Pt ( x ) = -L [ cos{I{t,x,rj))dr]. 
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Before proceeding further on, let us give a short description of the rest of the proof. We 
will estimate the integral (2.21) using the appropriate version of the Laplace method (see [Co65] 
for its description). In our case the application of the Laplace method meets some difficulties, 
since the expression under the integral contains two functions R and /. Therefore we introduce 
two intervals [—a, a] and [—/3,/3], on which R and I are controllable in terms of their Taylor's 
expansions. Then we split the integral into the sum of integrals over < a},{\i]\ G (a,/3]}, 
and {|?7| > /?}, and estimate these integrals separately. As in the standard Laplace method, the 
first two integrals are controlled by using Taylor expansion arguments. For the third integral, any 
standard considerations, like convexity arguments from [Fe77], cannot be applied. Therefore, we 
use the specific arguments based on the structure of the functional under consideration. 

Step 2: choosing a, j3. Following the explanations given above, we split the integral (2.21) into 
the sum 

1 / / I / ^ t ^+ i ^ x ))cos(I(t,x,r] + i^(t,x)))dr^ 



2tt 



\<a J\v\e(a,p] J\r)\>P J 



(2.22) 

= J ± (t,x) + J 2 (t,x) + J 3 (t,x), 
where a = a(t,x) and j3 = (3(t,x) are auxiliary functions. The function (3 is defined by 



(2.23) W ,x) = , "f^l 

V Mt(t,S(t,x)) 

Our aim in this step is to construct the function a in such a way that 

(2.24) < a(t, x) < Pit, x), (t, x) e A, 

Mz(t,S{t,x)) 



(2.25) 



M 2 (t,Z(t,x)) ^ a {t,x) ^ M*(t,Z(t,x)y 



a 3 (t,x) < 1 rr, 9(t) +x ->■ oo, (t,x)eA. 



M 3 (t,t(t,x)y 

By the Cauchy inequality and condition (Hi), we have 

M 2 3 (t, < M 2 (t, i)M,(t, « M 3 2 (t, 0, 9(t) + ^ oo, (t, G B. 
Hence, there exists a function k = K,(t,£), such that 

0, <t,t)<:M2(t,t)M: 1/2 (t,t), 

Without loss of generality, we can assume the function k to be locally bounded. Then we put 
(2.27) a(t, x) = cK(t, £(t, x))M 2 1/2 (t, £(t, x)) 

with some constant c > 0. By (2.26) and (2.11), we have (2.25). Since k is locally bounded, the 
constant c can be chosen small enough to provide (2.24). 

Step 3: estimating Ji(t,x) in (2.22). A straightforward computation shows that 

d d 3 d 2 

^R(t,X, V )\ v=0 = g^R(t,X, V )l =Q = 0, g^R(t,X, V )\ v=0 = -M 2 (t,t(t, X )), 
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(2.28) 
d 4 



R(t, x, 77) 



J J u 4 f 4 (t, s)e^ (t ' s)u cos (r]f(t, s)uj n{du)di 



<M 4 (t,£(t,x)), v eR, 



dr] 4 
which gives 

(2.29) - M 2 (t, £(t, x)) - ^M,(t, £(f , x)) < —R(t, x, r,) < -M 2 (t, £(t, x)) + ^M 4 (t, £(t, x)) 

for all i] G K. Therefore by the estimate for a 2 in (2.25) we get 
d 2 



(2.30) 



sup —R(t,x,ri) ~ -M 2 (t,£(t,x)), 

\ri\<a Of] 



inf — R(t,x,r}) M 2 (t,£,(t,x)), 0(t) + x ->■ 00, (t,x) e A. 



d 2 

\rf\<a d7] 

Next, similarly to (2.3) we get 
d 

I(t,x, V )\ v=0 = —I(t,x, V )\ v=0 



I(t,x,r))\ v=0 = 0, \—I(t,x,rj)\ < M 3 (t,£(t,x)). 



dr] 



Note that the equality for ^-1 holds true because z = i£(t,x) is a critical point for H{t,x,-). 
Hence the estimate for a 3 in (2.25) implies 

(2.31) sup |J(t,a;,77)| ->■ 0, 9(t) + x ->■ 00, (i,a;)eA 

r/|<ct 

Recall that if (t, x) = M 2 (t, f (t, x)) and D(t, x) = H(t, x, i£(t, x)) = R(t, x, 0). From (2.30) and 
(2.31) we get 



e R(t,x, v ) CQS 



\v\<a 



(2.32) 



I(t,x,r])dr]~ e R ^ 0) [ , 

271 e R(t,x,Q) f 

K(t,x) J lv 



dr] 



e 2 



\<y/K(t,x)a V27T 



dr) 



2tt 



e D{t ' x \ 0(t) + x^ 00, (t,x)eA, 



K(t,x) 

where in the last relation we used the lower estimate for a in (2.25). Thus, 



(2.33) 



Ji(t,x) 



e D{t ' x \ 9{i)+x ^00, (t,x)eA. 



x 



Step 4 : proving that J 2 (t,x) in (2.22) is negligible. On the set < f3}, the function R is 
controlled by its Taylor expansion. Hence for the integral J 2 (t, x) we can apply standard arguments 
of the Laplace method. 

By (2.29) we have for \r]\ < (3 

R(t, x, rj) < R(t, x, 0) - ^M 2 (t, C(t, x))r] 2 , 
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which, together with the lower estimate for a in (2.25), gives 

\J 2 (t,x)\< [ e R(t,*,v) dr] < e R(t,zfi) r e - M ^t xW d V 

K Z - 64: ) P D(t,x) 



e 2 drj <C Ji(t, x), 9(t) + x ->■ oo, (t,x)eA. 



\/K{t,x) J\y\> ay /K^x) 

Step 5: proving that Js(t,x) in (2.22) is negligible. By (2.3), 

\J 3 (t,x)\ < -L / e R ^d V 
^ J\ v \>p 

<J_ e D(t, x ) f exp(- / / e /( *' s) ^(l - cos(/(t, s)77u)V(du)dsl dry. 
2vr Jir,|>/3 I JiJr J 



Therefore, by (2.33), to prove J?,(t,x) <C Ji(t,x) we need to check that 



(2.35) / e- A ^ ri) dri<^K- 1/2 (t,x)=M2 1/2 (t,^t,x)), 6(t) + x ->■ oo, (t,i)eA 

•/|r?|>/3 

where 



A(t,x,77)= / [ e f ^ s) ^ x)u (l-cos(f(t,s)r]u))fi(du)ds. 
J i Jr 

Recall that £(t,x) > 0. Then for any a G (0, 1) we have 



A(t, x, 77) > / / «;/(*,«)«*,*)«(! _ cos(/(i, s)r]u))fi{du)ds 

(2.36) ii{(s,«)e/xi:/(t,s)u>o} 

> (1 - <7)9(t, 77, M + ) + <Te rx(t)€(t ' x) e(t, 77, [rx(0, +00)). 
Condition (H 3 ), combined with the trivial observation that Q(t,r],M, + ) is non-negative, yields 

(2.37) sup / e-( 1 - CT ) e ^ R+ )^<oo, 

teT Jr 

provided that a is chosen such that (1 — <r)(l + 5) > 1. Applying condition (H A ) with e = /3(t,x) 
gives 



e rx(*)«t,* 



0(t,77, [r X (t),+oo)) > ce r xW^9(t)(((3(t,x) X (t)) 2 A l), ^ > 0(t,x). 



Thus, in view of (2.37), to show (2.35) it is enough to prove 
(2.38) 

exp ^-ae rxmt ' x) 9(t) x) X {t)f Al)]« A^ 1/2 (t, f (*, x)), 0(t) + x ->• 00, (t, rr) G .A 

for every a > 0. By the definition (2.23) of f}(t,x), we have 

(2.39) ((/3(t,x) X (t)) 2 Al) = ((x- 2 (0^||) Vl)" 1 . 
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Condition (H 2 ) and the assumptions on 9 and Xi imposed prior to Theorem 2.1, yield for any 
a > 

(( X " 2(t) ^l|) vl ) ln - M 2(t,0«^K xW5 , 0(i) + £^oo, 

This relation, combined with (2.39), (2.11), and the relation £(t,x) > 0, yields 

In M 2 (t, £(t, x)) < ^(i^***^ x)x(t)) 2 A l) , 9(t) + x ->• oo, (t, x) G .A, 

which in turn implies (2.38) and completes the proof of (2.35). 
We have proved 

J 2 (t, x) < Jx(i, x), J 3 (t, x) < Ji(t, x). 
By (2.33) we get the statement of the theorem. □ 

3. Explicit conditions: fixed time setting 

Our further aim is to give explicit and tractable sufficient conditions which provide assumptions 
(Hi) - (H 4 ) of Theorem 2.1. In this section we consider the case where the time variable is fixed. 
Therefore everywhere below in this section we assume 

T = {t}, A = B = {t}xR, T={t}. 

We skip the variable t in the notation and write, for instance, Y, f(s), M.k(Q instead of Y t , f(t, s), 
Mk(t,£), respectively. 

In the fixed time setting the assumptions (Hi) - (H4) look more simple: in particular, functions 
6(t) and x(t) degenerate to some constants 9 and x- Therefore it is appropriate to introduce the 
set of conditions which will be useful later on. 

(Hi) M 4 (£)«.M|(0,£^+oo. 

(H 2 ) In V l) + \n\nM 2 (0 « ^ +00. 

(H 3 ) There exist R > 0, 5 > such that 
(3.1) e(z,R + ) > (l + 5)\n\z\, \z\>R. 

(H4) There exist q > and d > such that for every e > 

inf 0(z, k+00)) >^fe 2 AlV 

|z|>£ V / 

One can easily see that in the fixed time setting conditions (Hi) - (H4) are equivalent to (Hi) 
- (H4). Indeed, the constants 6 > and x > 0, which come, respectively, from the functions 
9(t),x(t), are suppressed in (H 2 ) by the term £. In (H4), the constant x can be eliminated by a 
proper change of the constants r and 9; we denote these new constants by q and 

Clearly, F is infinitely divisible with the Levy measure 



^f( A ) = l u f( s ) eA n(du)ds. 

J J IxR 
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In what follows we demonstrate that conditions (Hi) - (H4), which are in fact the assumptions 
on /if, can be verified efficiently in the terms of the kernel / and the initial Levy measure ji. 



3.1. Assumptions (Hi) and (H 2 ). Observe that (Hi) and (H 2 ) control the growth rate of the 
"tails" of /if. The following two lemmas show that these assumptions can be verified in the terms 
of similar "tail" conditions imposed on Denote 

(3.2) ^i(0 — I u(e^ u - l)n(du), M k (£) = [ u k e^ u fi(du), k>2. 

Jr Jr 

Clearly, 

(3.3) M k (i) = Jj k (s)M k (f(s)i)ds, k>\. 

Lemma 3.1. Assume 

(Ti) there exists 7 G (0, 1) such that M 4 (£) < M|(7f), £ ->■ +00. 
Then (Hi) holds true. 

Proof. Under the assumption (1.11) we have M k (£) —¥ +00, £ — > +00, k > 1. In addition, by 
Holder inequality, for every a G (0, 1) and k > 2, 

~l l-a 



M k (ai) < [M k (0] a 



u k jj,(du) 



implying 

(3.4) M fc «)<M fe (0, 

Denote F = esssup s6J /(s) and 7/ ;7 = {s : /(s) > 7T 1 }. Recall that / is assumed to be bounded, 
which together with (2.5) yields F G (0, +00). Since f 2 is integrable on I and / is bounded, 
f k , k > 3 is integrable as well. Then (3.4) yields 



(3.5) 



M k (0 



f k (s)M k (f(s)0ds) l + o(l) 



/,7 



, £^+00. 



Since M 2 is convex and M 2 (£) — > +00, as £ — > +00, there exists £ such that M 2 is increasing on 
[£o, +00). Then for £ > CoT -1 ^ -1 we have 

/ f 2 (s)M 2 (f(s)0ds\ = [ j f 2 (si)f 2 (s 2 )M 2 (f(s 1 )i)M 2 (f(s 2 )i)dsids 2 

>M 2 2 ( 7 F0^ f(s)ds 

Similarly, for sufficiently large £ we have 

f f\s)M 4 (f(s)i) ds < M 4 (F£) /" f\s) ds. 

These relations, together with (3.5) and (Ti), imply (Hi). □ 
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Lemma 3.2. Assume 

(T 2 ) In (m, Vl)+ lnlnM 2 (0 « ^ +00. 



Then (H 2 ) holds true. 
Proof. Fix an arbitrary 7 e (0, 1). For £ large enough, we have by (3.5) 

-M 2 (0 ~ / f(s)M 2 (f(s)0 ds < M 2 (F0 ( / / 2 (s) cfe 

which together with (T 2 ) gives 

(3.6) lnln.M 2 (0<£, £^+00. 

On the other hand, (T 2 ) implies that for every e > and £ large enough, 

M 4 (£) < e £ «M 2 (0- 
Then for £ large enough we have by (3.5) 

M 4 (0 ~ / / 4 (s)M 4 (/(s)0 < e £F « / / 4 ( S )M 2 (/(s)0 ds 
< FV F « / / 2 ( S )M 2 (/( S )0 rf S ~ *V F *M 2 (0- 



Consequently, 



Since e > is arbitrary, this relation combined with (3.6) implies (H 2 ). □ 

3.2. Assumptions (-£^3) and (H 4 ). It will be convenient to consider, together with the assump- 
tion (H 3 ), its stronger version 

(HI) Q{z,R + ) > ln|z|, z^oo. 

To proceed with the assumption (i^J), we introduce several conditions on the kernel /. 
(F 1 ) /,(/(*) V0) 2 rf S >0. 

(F 2 ) On some interval [a, 6] C /, the function / is positive and has a continuous non-zero 
derivative. 

(F 3 ) On some interval (—00, 6] c /, the function / is positive, convex, and has at most expo- 
nential decay at —00; that is, there exists 7 > such that 

(3.7) lim e" 7 7(s) = +00. 

(F 4 ) On some interval (—00, b) C /, the function / is positive, convex, and has a subexponential 
decay at —00; that is, (3.7) holds true for every 7 > 0. 

Note that when % increases from 1 to 4, the respective conditions (Fj) become stronger. Con- 
dition (Fi) is just our standing non- degeneracy assumption (2.5), listed here for further reference 
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convenience. Conditions (Fi) - (F 4 ) are well designed to handle the particularly interesting classes 
of processes, mentioned in the Introduction. Namely, 

(1) for the Levy process Z, one has fit, s) = Ir 0jt i(s), which satisfies (Fi) for every t > 0; 

(2) for the non-stationary version of a Levy driven Ornstein-Uhlenbeck process, one has 
f(t, s) = e 7 ^ _s ^l[o,t](s), which satisfies (F 2 ) for every t > 0; 

(3) for the stationary version of a Levy driven Ornstein-Uhlenbeck process, one has f(s) = 
e 7S l(_ OOi0 ](s), which satisfies (F 3 ) with 6 = 0; 

(4) for the fractional Levy motion, one has f(t,s) = fpj^ya) — s )+ ^ 2 — ( — s ) + ~^ 2 ] ' 
which for every t > satisfies (F 4 ) with 6 = 0. 

Recall several conditions which appeared in the literature in the context of the problem of 
studying local properties of infinitely divisible distributions. 

A Levy measure v on R is said to satisfy the Hartman-Wintner condition ([HW42]), if 



(3.8) / (1 — cos zu)v(du) >> In \z\, 

Jr 



z — > oo. 



Clearly, (i7J) is exactly the assumption on /x/, restricted to M + , to satisfy the Hartman-Wintner 
condition. 

An elementary inequality 

(3.9) cx 2 l\x\<i < 1 - cosx < x 2 A 1, i6l 

(where c > is some constant) provides the following pair of conditions, sufficient and necessary 
for (3.8), respectively: 

(3.10) / (uz) 2 is(du) > In \z\, \z\ ->■ oo; 

J\u\<\z\-i 

(3.11) / [{uz) 2 A l\v{du) 3> In \z\, \z\ ->• oo. 

Jr 

Condition (3.10) was introduced in [Ka81], and is called the Kallenberg condition. Condition 
(3.11) was introduced in [Ku06], where it was proved to be necessary for the existence of a bounded 
transition probability density of the solution to a (not necessarily linear) Levy driven SDE. At the 
same time, for an Ornstein-Uhlenbeck process (1.7) with non-trivial drift (7 7^ 0), this condition 
is sufficient for the existence of C°° distribution density ([BK08]). Thus, for the non- stationary 
version of the Ornstein-Uhlenbeck process (1.7), condition (3.11) is a criterion. 

We denote by //+ the restriction of \i to M + , and formulate the following set of "non-degeneracy" 
conditions on the measure /i. 

(Ni) /!+ satisfies (3.10). 

(N 2 ) /i+ satisfies (3.11). 

(N 3 ) pl(R + ) = +00. 

(N 4 ) fi(R+) > 0. 
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Note that when % increases from % — 1 to % — 4, the respective conditions (iVj) become more 
mild; (AT 4 ) is just our fixed non- degeneracy assumption (1.11), listed here for further reference 
convenience. 

Lemma 3.3. Assume for some i — 1, . . . ,4 conditions (JVj) and (Fj) /io/d. 
T/ien (H|j /io/ds true. 

Proof. Case i — 1. From the positivity of 1 — cosx and the first inequality in (3.9), it follows that 



9(^, M + ) = / / (1 - cos(uf(s)z))n(du)ds 

' (s,u):uf(s)>0 



> / / (1 — cos(uf(s)z))n(du)ds 

J J (s,u):u>0,0<uf(s)<l/\z\ 

>c u 2 f 2 (s)z 2 fj,(du)ds>cz 2 / f + {s) 2 ds I u 2 fi(du) : 

J J(s,u):u>0,0<uf(s)<l/\z\ [J I J ^(O^FIzl)- 1 ) 

here we keep the notation F = esssup sg/ /(s). Combined with (3.10) for the estimates above 
provide (Hf). 

Case i = 2. Since / is positive on [a, b], we have 

G(z, R + )> (1 - cos(uf(s)z)) n(du)ds. 

J J(a,b)xR+ 

Let us show that 

(3.12) / (1 — cos(x/(s))) ds > c(x 2 A 1) 

J a 

holds true with some constant c > 0, which would imply (Hf) provided that the assumption (3.11) 
is satisfied. Consider the function 



T(x) = j (1 - cos(x/(s))) ds. 

J a 



Clearly, T(x) ~ c\% 2 as x — > 0, with c\ = (1/2) f 2 (s) ds > 0. Further, one can write 

rf(b) 

(3.13) T(ar) = / (1 - cos(xv))g'(v) dv, 

where g := By our assumptions on / we have g G C 1 , which implies T(x) > for every 

x 7^ 0. Finally, by the Riemann-Lebesgue lemma, 

/•/(?») 

/ cos(xv)g'(v) dv — > 0, rr — >■ oo, 

which implies lim^oo T(x) > and completes the proof of (3.12). 
Cases i = 3 and i — 4. We show that the inequality 

(3.14) / (1 - cos(rr/(s))) ds > cln|x| 

J — oo 
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holds true (i) for some c > and \x\ large enough provided that / satisfies (F3); (ii) for every 
c > and \x\ large enough provided that / satisfies (F4). Keeping the notation g for the inverse 
function for /, we have 

fb rf(b) 
T(x) :— / (1 — cos(x/(s))) ds — / (1 — cos(xv))g'(v) dv. 

J -00 Jo 

Since / is convex, /' is non-decreasing. In addition, / itself is increasing: this follows from 
the convexity, condition (3.7), and the fact that f(s) — > as s — > —00 (which comes from the 
integrability of f 2 (s)). Therefore, g'{v) = [/'(fl , ( t '))]~ 1 i s positive and non-increasing. 

By positivity of g', 

/■/(b) 

(3.15) T(x) > / (1 — cos(xv))g'(v) dv 



J*/(2\x\) 

when 7r/(2|x|) < f(b). Denote Ik ■= ^ 2 , ,k >1. Then, since g' is positive, 

{-if I cos(xv) g'(v)dv > 
for every k > 1, and since g' is non- increasing, we have 

(3.16) / cos(xv)g'(v ) dv + / cos(xv)g'(v) dv < 

•J Ik-l J h 

for every even k > 2. Note that, on the axis [0, +00), the "negative" interval Ik-i is located to 
the left from the "positive" interval Ik- Then, for any A > 0, inequality (3.16) still holds true 
with and Ik replaced, respectively, by Ik-i H [0, A] and fl [0, A]. Consequently, for any 
A>tt/(2\x\) 

pA 00 „ 

/ cos(xt>)5f / (f ) dv = / cos(a;f )g'{v ) dv 
A/(2|x|) fe=1 i/ fc n[o,A] 

= I / cos(xf )(?'(f ) dv + cos(xv)g'(v) dv\ < 0. 

Therefore we obtain by (3.15) 

rf(b) 

(3.17) T(x) > / g'(v) dv = g(f(b)) ~ ^(tt/(2|x|)) = 6 - g(*/{2\x\)) 
for |x| large enough. It follows from (3.7) that 

g(v) 



p := liminf 



v^o \ \n(l/v) 

is positive when / satisfies (F 3 ), and equals to +00 when / satisfies (F4). Combined with (3.17), 
this yields (3.14). 
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Now we can complete the proof. In the case % = 3, take c > and Q > such that (3.14) holds 
true for \x\ > Q. Since /i(R + ) = +oo, there exists q > such that fJ,([q, +oo)) > (1 + 5)c~ 1 . Then 
(3.14) with x = uz implies 
(3.18) 

Q(z,R + )> ( (1 - cos(uf(s)z)) ds\ n(ds) > cp([q, +oo)) \n(\qz\), |^| > g" 1 ^, 

J[q,+oo) \J-oo J 

which provides (-fff) because ln(|gz|) ~ ln|z|, \z\ — > 00. 

In the case i = 4, the assumption /x(M + ) > implies the existence of q > for which 
/i([g, +00)) > 0. Take c satisfying c/i([g, +00)) > (1 + 5), and let Q > be such that (3.14) 
holds true with this c and \x\ > Q. Then (3.18) holds true as well, which provides (-fff). □ 



Lemma 3.3 shows that the kernel / is "smoothifying" in the following sense: when / satisfies 
some additional assumption like (F2) - (F4), the Hartman-Wintner type condition (H3) holds true 
under milder assumptions on the Levy measure of the noise. The following lemma shows that 
such "smoothifying" effect concerns the condition (#4), as well. 

Lemma 3.4. Under the assumption (1.11) assume additionally that the function f satisfies (F 2 ). 
Then (H4) holds true for q > small enough. 



Proof. Similarly to the proof of Lemma 3.3, case % = 2, we assume that / is positive on [a, b}. Take 
p > such that fi([p, +00)) > 0. Then, for < q < pmin se ( ajf) ) f(s), we have by (3.12) 

Q(z, [q, +00)) > / / (1 — cos(uf(s)z))dsfj,(du) 

J u>p J a 

(uz) 2 A l) n(du) > cp([p, +00)) [{pzf A l) , 



> c 

'u>p 

which implies the required estimate. □ 

To proceed with the assumption (H4) when / is not "smoothifying" , recall that a finite measure 
x is said to satisfy the Cramer's condition if 



(3.19) sup / e iyz x{dy) 

\z\>e JR 



< x(R) for all e > 



\z\>e 

(see, for example, [Lu79] or [IL65], chapter 3 §3). Cramer's condition means that x is in some 
sense regular. For instance, if x has a non-trivial absolutely continuous part, then (3.19) follows 
from the Riemann-Lebesgue lemma, although, in general, a measure satisfying Cramer's condition 
should not be necessarily absolutely continuous (see Example 3.3 below). 

Note that (3.19) leads to 

S(e) := sup / (1 — cosyz)x(dy) > for all e > 0. 

|z|>£ JR 
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In addition, assuming x to have finite second moment, we get S(e) ~ ce 2 as e — > with some 
positive c, and thus 

(3.20) 5(e) = sup / (1 - cosyz)x(dy) > c(e 2 A 1) for all e > 

|z|>£ </R 

and some positive c. Note that the function Q(z,A) involved in (i/4) is just the term under the 
supremum in (3.20), with x equal to pf restricted to A. By the standing assumptions on p and /, 
the measure pf restricted to R \ (—q, q) has finite second moment for any q > 0. Therefore, (#4) 
holds true, provided that for some q > the restriction of /i/ to [r, +00) satisfies the Cramer's 
condition. 

Lemma 3.5. Assume in addition to standing assumptions on p and f that 

(C) for some p > the restriction of p to [p, +00) satisfies the Cramer's condition. 
Then (H4) holds true for q > small enough. 

Proof. Take r < ^Fp with F = esssup se/ /(s) and some 7 e (0, 1) . Then 
@(z, [r, +00)) = // (1 — cos(uf(s)z))pi(du)ds 

J J (s,u):uf(s)>r 

> / (1 — cos(uf(s)z))/j,(du)ds > I / dsjEp^Fe), \z\>e, 

Jf{s)>jF Ju>p \Jf(s)>-yF J 

with "B p {e) = sup| z | >£ J p °°(l — cosuz)p(du). Since S p (e) satisfies (3.20) and the set {s : f(s) > 7F} 
has positive Lebesgue measure, we obtain the required estimate for <d(z, [q, +00)). □ 

To summarise, let us formulate in the fixed time setting the asymptotic results for the distribu- 
tion densities of particular processes, listed in the Introduction. 

Corollary 3.1. Let Y be a Levy driven stochastic integral, specified below. Assume that the Levy 
measure of the noise satisfies (1.11), (1.5), and "tail" conditions (T\), (T 2 ). 

Then for every t > the distribution density p t exists, belongs to C£° , and satisfies 
(3,1) PlW ._l_ e ^., 

with respective functions K(t,x) and D(t,x), in the following cases: 

(1) Y is the Levy process Z , p satisfies (N ± ) and (C); 

(2) Y is the non- stationary version of a Levy driven Ornstein-Uhlenbeck process, p satisfies 
(N2); 

(3) Y is the stationary version of a Levy driven Ornstein-Uhlenbeck process, p satisfies (N 3 ) 
(in that case, pt(x), K(t,x) and D(t,x) actually don't depend on t); 

(4) Y is the fractional Levy motion. 
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3.3. Examples. In this section we give several examples that illustrate the conditions on the 
measure /x, introduced above. 

The first two examples illustrate two typical situations where "tail" conditions (Ti) and (T 2 ) 
hold. 

Example 3.1. Let /J be supported in a bounded subset. Denote by a + the minimal positive 
constant a such that fi((cr, +00)) = 0. One can easily show that for all e > and k > 1 one has 

(3.22) M fe (0»e^+- £ * M fc (0 - ^MW})e CT+? « e CT+ «, £ +00. 

This relation yields both (Ti) and (T 2 ). Indeed, for £ > small enough one has 7 := 2 {a+-e) e -0 
and 

M 4 (0 < e ( CT++£ ^ = e 2 (-+- e )T5 < M|( 7 0, 
which is (Ti). Similarly, for £ large enough 

In Vlj+ InlnM 2 (0 < 2e£ + ln£, 

which provides (T 2 ) because £ > is arbitrary. 
Example 3.2. Assume that for u large enough 

(3.23) QR e ~^ " MK +00)) " Q ^ e ~ bU "> 
where b > 0, (3 > 1 are some constants, and Q is some polynomial. 

For a > denote 



J [<7, + Oo) 



Clearly, 

for any A > 0, and 

This means that, for any er > 



M fc (0 ~ / M fc e«>(rf M ) 
; 3.2 h J[a ' +Qo) 

= a h e^n([a, +00)) + / A;?/" 1 + £u k e ? >([w, +00)) du. 

./[<7, + 0o) 



For any a > 0, m G Z, we have 

(3.25) / M m e« u e- 6u ' 3 ~ Cl (/3, 6, m^^TOTV 2 ^^^ , £ ->■ +00, 

where c 2 (/3) = — (3~~^ (we have no need to specify the constant ci(/3, b,m)). One can prove 

(3.25) applying the Laplace method in a standard way; we omit the detailed calculations. 
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Take a large enough; then (3.23) holds true for u > a. Then (3.24) and (3.25) yield for every 
A- > 1 . 

(3.26) 1 e «»C0^A < Mk{0 < Q fc (£) e *(fl^€^ 

for £ large enough, where Qk is some polynomial. 
By (3.26), 

M 4 (0 < Q 4 (0e C2 ^^ « (-^f mb ^^) 2 < M ^ e +oo, 

g 

as soon as 27"- 1 > 1, which implies (7\). Further, for £ large enough (3.26) gives 

" 1 (ti vl ) +lnlnM2(0 

< In ((WOfttt) V l) + In (c 2 (P)b^ lnQ 2 (0) + ln£ « {, { -+ oo, 

which provides (T 2 ). 

The following example illustrates condition (C) and the relations between the conditions (N\) 
- (AT 4 ). All the measures in this example have bounded supports, therefore "tail" conditions (Ti) 
and (T 2 ) are satisfied. 

Example 3.3. (a) Let p = ^2™ =1 n p 5 n -i, P < 1; the assumption on p provides that // is a Levy 
measure. In the case p G (—1,1), the asymptotic behavior of the integrals f^ <£ u 2 fi(du) is the 
same as for the a-stable case with a — 1 + p, i.e. is of a power type: 

(3.27) / u 2 p(du) x e 2 ' a , e^O 

J\u\<e 

(cf. [068], [Pi96], [IK05]). Therefore, conditions (Ni) - (N4) holds true. The analogy with the 
a-stable case is not complete: condition (C) does not hold, because for every r > the restriction 
of p to [r, +00) has finite number of atoms. Therefore, statements (2) - (4) in Corollary 3.1 hold 
true, but one can not claim (3.21) for the Levy process Z itself. Statement (1) of Corollary 3.1 
becomes applicable when p is replaced by p + x, where x is a measure with a bounded support, 
satisfying Cramer's condition. For instance, either x may be absolutely continuous (and then 
Cramer's condition is provided by the Riemann-Lebesgue lemma), or x may be equal to the 
Cantor measure on [0, 1] (and then Cramer's condition is verified by straightforward calculations). 

When p = -1, (JVi) and (N 2 ) fail, but (N 3 ) and (jV 4 ) hold true. When p < -1, only (jV 4 ) hold 
true, while (Ni) - (jV 3 ) fail. It is clear that, in the latter case, the laws of the Levy process Z 
and of the non-stationary version of a Levy driven Ornstein-Uhlenbeck process contain non-trivial 
discrete components. Therefore one definitely can not expect any asymptotic relation like (3.21) 
to hold for these processes. On the other hand, (3.21) holds true for the fractional Levy motion 
Z H with H E (1/2, 1). This well illustrates the "smoothifying" role of the kernel /. 

(b) v = J2^=i ^(n!)- 1 - Then condition (iVi) fails, while (N 2 ) - (iV 4 ) hold true, see [Ku06], 
Example 2.3 or [BK08], Example 1. In these examples, it is shown that the law of Z t is singular 
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for all t > 0. Thus, the asymptotic relation (3.21) clearly can not be valid for the Levy process 
Z itself. In this case, the Levy measure provides some "hidden smoothness" in the sense that the 
law of the Levy process Z is singular, but the distributions of the respective (both non-stationary 
and stationary) Levy driven Ornstein-Uhlenbeck processes and fractional Levy motion possess C£° 
distribution densities which, moreover, admit asymptotical description (3.21). 

In the last example, in the case of a Levy process, we compare our conditions with those 
introduced in [CKK08] and [CKK10]. 

Example 3.4. In the paper [CKK08] the authors give the transition density estimates for a 
symmetric a-stable-like process whose jump intensity kernel J(x, y) is of the form 

where c(x,y) is a symmetric Borel measurable function on MP x R n , bounded from above and 
below by two positive constants. When c(x, y) = c(x — y), that is, J(x, y) = J(x — y), this process 
is a Levy one with the Levy measure /i(dx) = J(x)dx. We check in the one- dimensional case that 
such a Levy measure satisfies the conditions imposed above. 

Since the Levy measure \i has bounded support, the exponential integrability condition is sat- 
isfied and, moreover, conditions (7\) and (T 2 ) hold true (see Example 3.1). By the Riemann- 
Lebesgue lemma, the absolute continuity of ji implies condition (C). Finally, (3.27) holds true, 
which provides the Kallenberg condition (3.10) for the measure ji. Since /i is assumed to be 
symmetric, this yields (Ni). 

The paper [CKK10] is devoted to the estimates of the transition density of a Markov process 
whose jump intensity J(x,y) satisfies 

(3-28) — < J(x, y) < -. — rr, x,yeR n xR n , x^ y, 

\x — y\ n (p{C2\x — y\) \x — y\ n (p(c4\x — y\) 

for some q, % — 1,2, 3, 4, where : [0, oo) — > [0, oo) is of the form <f)(r) = 0i(r)'0(r), r > 0, and 

i) ip is increasing on [0, oo), ip{j-) — 1 for < r < 1, and for some < 71 < 72, f3 > 0, 

(3.29) Cl e 7ir ^ < V(r) < c 2 e< 2r \ 1< r < 00; 

ii) (f>i is strictly increasing on [0, 00) with 0i(O) = 0, 0i(l) = 1, and, in particular, satisfies for 
c 2 > ci > 0, c 3 > 0, < P 1 < f3 2 < 2, the inequality 

3.30 ci - < < c 2 - for every < r < R < 00. 

\r J 0i(r) V r / 

Again, let n — 1 and J(x, y) = J(x—y), where J is the density of the Levy measure /x. To achieve 
the exponential integrability (1.5) we need to assume /3 > 1 in (3.29). Assuming additionally that 
7i = 72, one has (Ti) and (T 2 ) (see Example 3.2). By (3.30), the Levy measure \x satisfies the 
lower bound in (3.27) with a = j3i; that is, 

u 2 fi(du) > ce 2 -? 1 

\U\<6 
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with some positive c and e > small enough, which implies (iVi). Finally, condition (C) holds 
true by the absolute continuity of fi. 

Since (Ni), (C), (Ti), and (T 2 ) hold true, by statement (1) in Corollary 4.1 and Corollary 5.1 
below, the transition probability density of the Levy process Z satisfies (1.12) and either (1.17) (in 
the "truncated" case [CKK08]) or (1.18) (in the case treated in [CKK10]). Let us compare these 
relations with the estimates for the transition probability density of a symmetric jump process 
from [CKK08] and [CKK10]. 

For t > t , these estimates are given in the form 
(3.31) C l9t (C 2 \x - y\) < p(t, x, y) < C 3 g t (C 4 \x - y\), 

where Ci, . . . , C4 are some positive constants, and 



-\x\\n s — j V ( ^exp 

with 5 = 1 in the "truncated" case [CKK08], and 5 = -^-^ in the case treated in [CKK10] id is 
the dimension of the space; in the current paper d = 1). 

For a Levy process, (3.31) with p(t,x,y) = p t (y — x) is closely comparable with (1.12) and 
(1.17), (1.18). When ^ is large, (1.17), (1.18) directly provide (3.31) with g t replaced by 



exp ( —\x\ In 



\x\ 

T 

On the other hand, one can show easily that on every bounded set the function JCz is bounded 
and bounded away from 0, and the function V z satisfies 

-dxx 2 < V z (x) < -d 2 x 2 

with positive constants d±, d 2 . Thus, when is bounded, (1.12) provides (3.31) with g t replaced 
I'.v 



1 f \x 



t^ exp y t 

Note that (1.17) and (1.18) are somewhat more precise than (3.31): by choosing | large enough, 
one can make the constants Ci,c 2 therein to be arbitrarily close to a given constant c*, while in 
(3.31) respective constants C 2 and C 4 are different and fixed. 

Although having a non-trivial intersection, the classes of processes, treated in our case and in 
[CKK08] and [CKK10], are substantially different. Our approach, based on the Fourier transform 
technique, is not applicable to the class of symmetric jump processes from [CKK08] and [CKK10] 
in the whole generality. On the other hand, this approach is applicable to particularly interest- 
ing processes which can not be studied by the technique of [CKK08], [CKK10], including non- 
symmetric Markov processes (like the Levy driven Ornstein-Uhlenbeck process) and non-Markov 
processes (like the fractional Levy motion). 



12 
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4. Explicit conditions: time-dependent setting 



Our further aim is to consider conditions of Theorem 2.1 in the general, i.e. time-dependent, 
setting. To make the exposition reasonably short, we address this problem in a particular case of 
the self-similar kernel /; that is, we assume that 



(4.1) 



f(t,s) = X (t)f 



e(t)J' 



teT, s el 



with some functions / : R — > R and x, 9 : T — > (0,+oo). Assumption (4.1) is satisfied for 
particularly interesting processes like the Levy process Z and the fractional Levy motion Z H . In 
these cases we have, respectively, 



(4.2) 



(4-3) /(*) 



f(s) = l m {s), X (t) = l, 0(f) =f; 



1 



r(# + i/2) 



H-l/2 



s) 



H-l/2 



X (t)=t H -^, 9{t)=t. 



For the function f(s) we keep our standard standing assumptions: it is bounded and satisfies 
(2.3), (2.5). For the Levy measure fx we assume (1.11) and (1.5) to hold true, as before. 

Denote, similarly to Section 2, 

e(z,A) = 



{(s,u)GRxR: f(s)uGA} 



(1 — cos(f(s)zu))[/,(du)ds, zeR, 



V(z) = [ [ (e~ izf ^ u - 1 + izf(s)u) fx(du)ds, z e C, 
Jr Jr. 



H(y,z) = iyz + y(z), M k {Q 



dC 



;*«), k>l, 1,61, Ce 



Denote by ((y) G R the unique solution to the equation 

d 



(4.4) 
and put 
(4.5) 



H(y,iC) = 0, 



V(y) = H(y, i((y)), K{y) = M 2 (((y)) = — 2 H{y, i() 



C=C(w) 



Denote r{t) = x(t)9(t)- Further in this section we assume 9 and x to be bounded on every 
segment [a, b] C (0, +oo), and to be bounded away from on the whole T. Clearly, the functions 
9,x i n (4-2) and in (4.3) with H > 1/2 satisfy these assumptions. In addition, we assume that 

(4.6) 9(t) ->+oo, In ((lnx(0) V l) < \n9(t), t ->■ +oo; 

in the cases (4.2) and (4.3) this assumption holds true. 
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Theorem 4.1. Assume that the measure \x satisfies (T±) and (T2). Assume also that \x satisfies 
one of the conditions (Ni) and, respectively, f satisfies one of the assumptions (F^), i = 1, ... ,4. 
In the case % — 1, assume additionally that fi satisfies condition (C ). 

Then for every t > the law of Y t has a distribution density p t G C£° , and for every t > 



(47) p ' W ~7wV^^ eW ~' t + ^°°' (^)^'».+-)x» + 

Remark 4.1. The expression on the the right hand side of (4.7) is self-similar in the sense that the 
variable x, rescaled by r(t), is involved in this expression only as an argument of given functions 
/C and T>. Note that the Levy measure fi is not assumed to have a self-similarity property, and 
therefore, in general, the family of distributions Y t , t > is not self-similar. Thus, although the 
assumption (4.1) on the kernel / itself does not provide self-similarity for the distribution densities 
of Y t , t > 0, it is powerful enough to provide self-similarity for the asymptotic relation for these 
densities. 

Proof. The relations below follow easily from the self-similarity assumption (4.1): 

(4.8) H(t, x, z) = 6(t)H (JL, X (t)z^j , M k (t, = X k (t)6(t)M k ( X (t)0, k > 1. 

By the first relation in (4.8), we can rewrite the relation (2.10), which determines £ = as 

This means that s °l ves (4-4) with y = x/r(t), and therefore 

x 



Z(t,x)= X - L (t)( 



r(t) 



Combined with (4.8), this relation gives 



D(t,x) = 6(t)V[^-), K(t,x)= x 2 (t)6(t)ic(^-)=^1C / ' 



At) J' V ' ^ * w w \r(t)J 6(t) \r(t). 

Thus (4.7) would follow from (2.13) with A = [t , +00) x M + , provided that conditions (Hi) - 
(H4) are verified. 

In Section 3 we proved that under assumptions imposed on the Levy measure \i and the function 
f(s), conditions (Hi), (H 2 ), (Hf), and (H 4 ) hold true. Now we show that these conditions yield 
(Hi) - (H 4 ) with T = [t , +00), B = [t , +00) x R+, and with the function 6(t) replaced by $6(t) 
(the constant d comes from (H 4 )). 

The second relation in (4.8) gives 

uq] M,(t,Q 1 M 4 (x(t)Q 

{ -> Ml(t,0 9(t)Ml( X (t)0' 

Observe that, under our assumptions on 6 and Xi 

(4.10) t + £->oo implies 9(t) ->■ +00 or -» +00. 
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Therefore, (#1) follows from (Hj and (4.9). 
By the second relation in (4.8), 

which together with (H 2 ) and (4.10) gives 

ln ((^ 2(t) ^f) Vl ) « ln ^)+^> t + ^+oo. 

Similarly, 

In ((lnM 2 (t,0) V l) = In (in (x 2 (t)6(t)M 2 ( X (t)0) V l) 

= In ((ln X 2 (t) + \n9(t) + \nM 2 ( X (t)0) V l) ■ 

By (H 2 ), (4.10) and (4.6) one has 

In ( ( In M 2 (t, 0) Vl)« In 6(t) + t + ^+oo. 

This completes the proof of (H 2 ). 

By (i/f), for every x > there exists Q > such that 

Q{z,R + ) > x\n\z\, \z\ > Q. 

By the self-similarity assumption (4.1), we have 

e(t,z,A) = e(t)e( x (t)z,-$-rA). 



x(ty 

Denote 6* = inf 4 8(t), x* — m U x(t)- Then taking x = 6~ l {l + 5) and R = x* X Qi we obtain (#3). 
Finally, by (H4) we have 

inf Q(t,z,[qx(t) < +00)) =0(t) inf Q( X (t)z, [q < +00)) 
teT, |z|>e tsT, |z|>£ 



= 9{t) inf 9(z', [q < +00)) > W(f ) (( X (t)^) 2 A l) . 

t&T, z' >y f £ V / 



ter,| 2 '|>xW£ 

Thus, (#4) holds true with r = q and #(t) replaced by $6(t). Clearly, such a change of the function 
6{t) does not spoil conditions (Hi) - (iJ 3 ) proved above. □ 

Corollary 4.1. Assume the Levy measure of the noise satisfy (1.11), (1.5), and "tail" conditions 
(T l ), (T 2 ). Then 

(1) For the Levy process Z, assuming additionally \i to satisfy (Ni) and (C), one has (1.12). 

(2) For the fractional Levy motion Z H , one has (1.13). 
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5. Explicit asymptotic expressions as x ->■ +00 

Theorem 2.1, Corollary 3.1, and Theorem 4.1 describe the asymptotic behaviour of a distribution 
density precisely, but in an implicit form: functions K(t, x), D(t, x), JC(x), T>(x), involved in (2.13), 
(3.21) and (4.7), are defined in terms of the solutions to equations (2.10) or (4.4). In this section 
we study the asymptotic behavior of these functions as x — > +00, and deduce explicit asymptotic 
expressions for the distribution densities. 

In what follows, we mainly discuss the behavior of the functions K,(x) and T>(x) under additional 
assumptions on the Levy measure \i\ without any essential change of the argument, similar results 
can be obtained for the functions K(t, x), D(t, x) with a fixed variable t. To simplify the argument, 
we assume in the sequel / to be non-negative. This assumption is satisfied, for instance, for the 
Levy process and the fractional Levy motion (respective functions / are given in (4.2) and (4.3)). 
To shorten the exposition, we restrict ourselves to the cases where the Levy measure \i is either 
"truncated" (i.e. supported in a bounded set, see Example 3.1) or "exponentially damped" (i.e. 
having its "tails" satisfying (3.23), see Example 3.2). 

We keep the notation introduced in Example 3.1, Example 3.2, and Section 4; in particular, 
F = essup f(s), and <x+ is the extreme right point of the support of ji. 

Theorem 5.1. Assume the kernel f(t,s) to be of the form (4-1) with 9 and x satisfying (4-6). 
Assume that the measure fi satisfies one of the conditions (N^), and the respective function f(s) 
in (4-1) satisfies one of the assumptions (Fi), i — 1, . . . ,4. In the case % — 1, assume additionally 
fi to satisfy condition (C). 

1. If ji is truncated, then for any constants C\ > l/(a + F) and c 2 < l/(a + F) there exists 
y = y(ci, c 2 ) such that, for x/r(t) > y, 

(51) h (m) ln (4)) £ Mx) £ exp h (m) ln (m)) • 

2. If /i is exponentially damped, then for any constants 

13-1 

c 2 < (/r^ _ p-T^j b^F" 1 < ci 

there exists y = y(ci, c 2 ) such that, for x/r(t) > y, 

(5.2) exp (- Cl (£) .„¥ (^)) <p (W < exp (-e 2 (£) .„¥ (^)) . 

Proof. We consider in detail the case of a truncated Levy measure, and then outline the changes 
in the proof that should be made in the case of an exponentially damped Levy measure. 

Denote 

Mo(C) = *«) = / / (e C/(s) " - 1 - (f(s)u) ^{du)ds. 
Jr Jr 

Similarly to (3.3), one has 

M)(C)= / M (f(s)()ds, M (0:= / (e^-l-^u)fi(du). 
Jr Jr 
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To describe the asymptotic behavior of /C, T>, we need to analyze the behavior of Ai k , k — 0, 1, 2. 
For this, we analyze first the behavior of M k , k — 0, 1, 2. 

One can easily see that (3.22) holds true for /c = as well. From (3.22) we have for any k > 

(5.3) M*(0 ~ o*M (0, C^+oo. 
Moreover, the first relation in (3.22) provides that for every e > 

(5.4) M k (()~ [ f k (s)M k (f(s)()ds, C^+oo, 

Jf(s)>F-e 

(recall that we assume / to be non-negative), which together with (5.3) yields 

(5.5) M k (() ~ F k a k + M ((), C^+oo. 

Recall that K{y) = M 2 (((y)), and 

v(y) = -y((y) + Mo(((y)). 

The function ((y) is defined by the equation A4i(((y)) = y and, under our fixed assumption 
(1.11), we have ((y) — > +oo as y — > +oo. Hence, by (5.5), 

(5.6) K{y) ~ Fa + y, M (((y)) ~ (l/Fo+)y, y -> +oo. 
The second relation in the above formula yields 

(5.7) ~ -yC(y), y -> +oo. 

Similarly to (5.5), one can deduce from (3.22) that for any e > 
and consequently 

(5.8) ((y) ~ -J—\ny, y ->• +oo. 

Let us prove the lower bound in (5.1), the proof of the upper bound is similar and omitted. It 
follows from (5.7) and (5.8) that for any c > l/(a + F) we have for x/r(t) large enough 



(recall that r(t) = 6(t)x(t)). 

Since // is supported in a bounded set, it satisfies "tail" conditions (Ti), (T 2 ) (see Example 3.1), 
and thus we can apply Theorem 4.1. By Theorem 4.1 and (5.9), to prove the first inequality in 

(5.1) it is enough to take c G (l/(a + F), c±j and prove that for x/r(t) large enough, 



(5.10) —J Jri >expf(c- Cl )f4TV n ^ X 



r(t)\2KlC(x/T(t))- "V" J \x(t)J \r(t) 
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By (5.6) and (4.6), we have for any q > 0, 



frnn 1 ./ e JJ} > 1 P -g»/T(t) 

1 J r(t)]j 2nlC(x/T(t)) ~ 0(t) 

for x/r(t) large enough. On the other hand, for a fixed y > 1 and x/r(t) > y, 

exp ((c - ci) In (^y)) < e-^W^W) with g = ( Cl - c 2 ) my > 0. 

Hence (5.10) follows from the inequality 

lg-?6 _ e -lna-g6 y e ~q(a+b) > e ~1 ab 

a 

valid for a, 6 large enough. 

Let us discuss briefly the changes that should be made when the measure \x satisfies (3.23). 
Clearly, for any a > and k > j, 

MZ(0>a k - j M?(Z), £>0 
(see the notation in Example 3.2). This means that instead of (5.3) and (5.5) we have now 

(5.12) M fc (0»M,-(0, e^+oo, 

(5.13) MkiO^MjiC), C^+oo 

for any k > j. The latter relation with k — \,j — yields (5.7). From (3.26) and (5.4) it follows 
that for every e > for ( large enough 

e (C-e)cA < Ml (C) < e^ +e ^ , 

/ i_ g_\ _j_ g 

where = 1/3 z 9 - 1 - /3 ^ b^Ff-^. Consequently, 



(5.14) (( y )^^ — \nyj , V(y) ~ -y ^ — \nyj , y -> +oo, 

which means that the analogue of (5.2), with e d ^ v ^ x ^ T ^ instead of p t (x), holds true, and the only 
thing we need to verify is that the term 



6(t) 



r(f) y 2nlC(x/T(t)) 
is negligible. Note that this term is bounded: 



1 6(t) 1 
sup —-:\ -t. — = sup — = < +oo 

t l r(t) V 2tt t y ^2-K X 2 {t)6{t) 

because 6 and x are assumed to be separated from 0, and by (5.13) 

JC(x/r(t)) = M 2 (((x/r(t))) » Mi(ax/r(t))) = x/r(t), x/r(t) +oo. 
This observation provides the upper bound in (5.2). 
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On the other hand, it follows from (3.26) that 

In — ^ < £, £ ^ +00 
(cf. (5.12)). Similarly to the proof of Lemma 3.2, one can deduce from this relation that 

and consequently 

In /C(y) < C(y) + In y, y ->• +00. 
Together with (4.6) and the first relation in (5.14), this implies (5.11). Repeating the argument 
after (5.11), we obtain the lower bound in (5.2). □ 

For the Levy process Z and the fractional Levy motion Z H , Theorem 5.1 gives the following. 
Denote 

(5.15) c* = 1/(7+ 
in the case of the truncated Levy measure //, and 

g-i 

(5.16) c*=(/T^-/T^) b~~? 
in the case of the exponentially damped Levy measure [i. 

Corollary 5.1. Assume that the Levy measure satisfies (N\) and (C), then for the distribution 
density of the Levy process Z the following estimates hold. 

1. If fi is truncated, then for any constants c\ > c* and C2 < c* there exists y = y{c\,C2), such 
that for x/t > y (1.17) holds true. 

2. If fi is exponentially damped, then for any constants c\ > c* and c 2 < c* there exists 
y = y(c 1 ,c 2 ), such that for x/t > y, (1.18) holds true. 

Corollary 5.2. Assume that the Levy measure satisfies (1.11), then for the distribution density 
of the fractional Levy motion Z H the following estimates hold. 

1. If fi is truncated, then for any constants c\ > c* and Ci < c* there exists y = y(c\,C2), such 
that for x/t H+1 / 2 > y (1.19) holds true. 

2. If fi is exponentially damped, then for any constants C\ > c* and c 2 < c* there exists 
y = y(ci,C2), such that for x/t H+1 ^ 2 > y (1.20) holds true. 

We have mentioned in the beginning of the section that for fixed t the functions K(t, x), D(t, x) 
can be analyzed in the same way as K,{x) and V(x). Respectively, the analogue of Theorem 5.1 
can be proved for the density Pt(x) with fixed t without the self-similarity assumption (4.1). Let 
us formulate one statement of such a kind. 

Consider the stationary version X of a Levy driven Ornstein-Uhlenbeck process, and assume that 
\i satisfies (N 3 ). Then the distribution of X t , in fact, does not depend on t, and by Proposition 2.1 
has a C°° distribution density p, which we call the invariant distribution density of the respective 
Ornstein-Uhlenbeck process. Moreover, assuming the "tail" conditions (Ti) and (T 2 ) to hold, we 
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have by Corollary 3.1 the asymptotic relation for this density, which after trivial transformations 
can be written in the form 

(5.17) p(x) ; 1 e v{x \ x -> +oo, 

where /C, T> are defined by (4.5) with f(s) = e 7S l s <o- Similarly to Theorem 5.1, one can deduce 
from (5.17) the following statement (the proof is omitted). 

Proposition 5.1. Assume that the Levy measure fi satisfies (1.5) and (N 3 ). Then for the invari- 
ant distribution density of the Ornstein-Uhlenbeck process (1.7) the following estimates hold. 

1. If fi is truncated, then for any constants C\ > c* and c 2 < c* there exists y = y(ci,c 2 ), such 
that for x > y, 

exp (—cix In x) < p(x) < exp (— c 2 x In x) . 

2. If \i is exponentially damped, then for any constants c\ > c* and c 2 < c* there exists 
y = y(ci, c 2 ) such that for x > y, 

exp c\xhi~ x^j < p(x) < exp ^— c 2 x\n~ x^j . 
Here c* depends on /x only, and is defined, respectively, in (5.15) or (5.16). 

As we mentioned in the Introduction, there is a particular theoretical interest in studying the 
ratio (1.14) of the values of the invariant distribution density. One can see that the statement of 
Proposition 5.1 is not strong enough to provide an exact estimate for the ratio (1.14) because of 
different constants c\ and c 2 , involved in respective estimates. In the theorem below we provide 
the exact estimate for the ratio (1.14). 

Theorem 5.2. Assume that the Levy measure ji satisfies (T\), (T 2 ), and (N 3 ). 
Then for every bounded set A cR 

(5.18) r a (x) ~ e~ aax \ x -> +00, 
uniformly in a £ A. 

In particular, for any constants c\ > c* and c 2 < c* there exists y = y(ci,c 2 ,A), such that for 
x > y, a G A, 

(5.19) x- cia < r a (x) < x-° 2a 
when ji is truncated, and 

(5.20) x- cialn ^ x < r a (x) < x~ C2aXn ' T x 

when jj, is exponentially damped. Here c* is defined respectively in (5.15) or (5.16). 



Proof. By the inverse function theorem, 



(5.2!) A< W 



d 



Tc M l{0 



= \M 2 (((x)]-K 



34 



Then 



rfx 1 } M 2 (((x))dx^ ) Ml({(x)Y 



If /i is supported in a bounded set, then J\4s(C) ~ (l/cr + )A^4(C), C ~~ * +°° ( see (5-5)). If /i is not 
supported in a bounded set, then A^C) <C A^C), C — > +°° ( see (5.13)). In both cases, we have 

— In K,(x) — > 0, rr — >■ oo 

ax 

because (7\) provides (#1) (Lemma 3.1). Thus, for any bounded set A, 
K{x + a) ( f x+a d 



( f x+a d \ 
ex P I / ~r m ^-(y) dy I — > 0, rr — >■ +oo 



/C(ar) 

uniformly in a G A Therefore, by (5.17), 

uniformly in a G A. 

We have 



V(x + a) - = -a((x) - (x + a)[C(x + a) - C(s)] + -Mo(C(& + a)) - M (((x)). 
d_ 



Since £M (t) = -Mi(C) and M^Ciy)) = y, we get 



fx+a px+a 

(x + a) [((x + a) - C(x)} -M (((x + a))- M (((x)) = (x + a) I ('(y) dy - / y('(y) dy 



X 



px+a px+a px+a pa pa 

= / (x + a - y)C(y) dy = / / ('(y)dvdy= / / ('(x + r)dsdr. 

J x J x J y Jo J r 

Since \x satisfies (N 3 ) we have (1.11), and therefore A1 2 (C) — > +°°> C — > +°°- By (5.21), this 
yields 

C(x) -» 0, x ->■ +oo. 
From the above relations we deduce that 

V(x + a) — T>(x) —a((x), , x — > +oo 

uniformly in a G A, which completes the proof of (5.18). 

From (5.18) and (5.8) we deduce (5.19). From (5.18) and the first relation in (5.14) we obtain 
(5.20). □ 
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